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Abstract. We investigate the idea that the effect of the truncation applied in the TCSA method 
on the spectrum coincides with the effect of a suitable changing of the coefficients of the terms 
in the Hamiltonian operator. The investigation is done in the case of the critical Ising model 
on a strip with an external magnetic field on one of the boundaries. A detailed quantum field 
theoretical description of this model is also given, and we propose a description as a perturbation 
of the infinite coupling limit. The investigation is also carried out for a truncation method which 
preserves the solvability of the model. The results of perturbative and numerical calculations pre- 
sented support the above idea and show that the qualitative behaviour of the truncated spectrum 
as a function of the coupling constant depends on the truncation method. 
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1 Introduction 

The main subject of this paper is an investigation of the TCSA (truncated conformal space 
approach) method, which is a numerical method introduced in [T] for the calculation of the 
spectra and eigenvectors of Hamiltonian operators of perturbed conformal field theories in 
finite volume in two spacetime dimensions. The Hamiltonian operators of these theories 
have the form H = Hq + Hj, where Hq has a known discrete spectrum. The method consists 
in (numerically) diagonalizing the finite matrix H'^^^^ijic) = P {ric) {Hq + H j) P [ric) , where 
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P{nc) is the projection on the subspace spanned by the eigenstates of Hq belonging to the 
lowest Uc+l energy levels, and Uc is the truncation parameter. This gives an approximation 
to the eigenvalues and eigenvectors of H. 

The uses of the data obtained by TCSA include the verification of results obtained by 
other methods, for example in |2j, the extraction of resonance widths [3|, the mapping 
of the phase structure of certain quantum field theories as in [HE], and the finding of 
renormalization group flow fixed points as in [6]. Our investigation is related to the use of 
TCSA for the study of boundary renormalization group flows (i.e. of flows generated by 
boundary perturbations) between minimal boundary conformal field theories defined on 
the strip [0, L] x M. A brief review of the areas where boundary conformal field theory and 
flows between them play important role can be found in [7j. 

By boundary flow we mean a one-parameter family of models, the parameter being 
the width (or 'volume') L of the strip. In the simplest case the parameter can be taken 
to be a coupling constant h instead of L and the models have the Hamiltonian operators 
H = Hq + hHj. Hq is the Hamiltonian operator of a boundary conformal held theory, h is 
allowed to vary from to oo or from to — oo. if/ is a relevant boundary held taken at a 
certain initial time. 

A particular problem of interest is that of finding values of h other than 0, called 
fixed points, where the model corresponding to Hq + hHj is a conformal field theory, and 
identifying these conformal field theories. The TCSA can be used for this purpose if the 
Hilbert space of the conformal field theory corresponding to /i = consists of finitely 
many (or countably many) irreducible representations of the Virasoro algebra. Boundary 
conformal minimal models satisfy this condition. We refer the reader to [HI [9l [IHl [HI [6] for 
boundary flows and the application of the TCSA to study them. 

Regarding the spectra at nonzero values of h, conformal symmetry can be recognized 
by the equal distances between neighbouring energy levels; the representation content can 
be identified by the degeneracy of the energy levels. It should be noted that perturbation 
theory and other methods can also be used (see e.g. [7l[l3l[lll[l5],[l6l[I7j) to explore flows. 
An important problem regarding the TCSA is that it gives approximate data, and our 
knowledge of the precise relation between this data and the exact spectrum is still limited 
(see [18j for already existing results). A good understanding of the effect of the truncation 
could be used to improve the TCSA data and to explain the qualitative features of TCSA 
pictures of flows, and generally it would make the results obtained using TCSA data better 
founded. 

An idea proposed by G.M.T. Watts and K. Graham [lH [20] recently is that the effect 
of the truncation on the spectrum is equivalent to a suitable change of the coefficients of 
the terms in if, i.e. the spectrum of H^'^'^^{h,nc) is equal, in an approximation at least, 
to the spectrum of 

H''{h,nc) = So{h,nc)HQ + Si{h,nc)Hj + S2{h,nc)Hi^2 + ■ ■ ■ , (1) 
which we shall call renormalized Hamiltonian operator. sq,si,... are suitable functions 
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and Hj 2^ ■ ■ ■ are suitable operators. Hj 2, ■ ■ ■ should be primary or descendant bulk or 
boundary operators. Our main purpose in this paper is to investigate the validity of this 
picture. 

We do this investigation in the case of a relatively simple model: the perturbed bound- 
ary conformal field theory on the strip [0, L] x R with Hamiltonian operator 

H=^Lo + hL~'/^^y2{x = L,t = 0) . (2) 

The unperturbed {h = 0) model is the c = 1/2 unitary conformal minimal model, i.e. the 
continuum limit of the critical Ising model, with the Cardy boundary condition on the 
left and 1/16 on the right. jLq is the Hamiltonian operator of this model. Lq is the 'zero 
index' Virasoro generator. L will be kept fixed at the value L = 1. The Hilbert space of the 
unperturbed model is the single c = 1/2, /i = 1/16 irreducible highest weight representation 
of the Virasoro algebra. (Here and later on the coupling constant and the highest weight 
are both denoted by h, but it should be clear from the context which one is meant.) The 
field 01/2 is the weight 1/2 boundary primary field on the right boundary, which is 
also known in the literature as the boundary spin operator [2H [22l [23]. The coupling 
constant h can also be regarded as a constant external boundary magnetic field, which is 
coupled to the boundary spin operator. The model ^ is also referred to as the critical 
Ising model on the strip with boundary magnetic field. The perturbation h(f)i/2{L,t = 0) 
violates the conformal symmetry, which is nevertheless restored in the h — > ±oo limit. It 
is known that in the /i — > oo limit the c = 1/2, h = 1/2 representation is realized; in the 
h — oo limit the c = 1/2, h = representation is realized (see e.g. [I3l [HI [23l [24]). This 
can be written in a shorthand form as 

(1/2, 1/16) + 01/2 ^ (1/2,1/2) (3) 
(1/2, 1/16) - 01/2 ^ (1/2,0) . (4) 

We have chosen the model ([2]) because it is integrable, furthermore it is relatively easy 
to handle; in particular the spectrum can be calculated analytically in terms of simple 
functions and the application of Rayleigh-Schrodinger perturbation theory is also relatively 
easy. Investigations in other perturbed conformal minimal models, especially in the case 
of the tricritical Ising model and generally in the case of a perturbation by the field 0{i3), 
are carried out in [25] (see also [20]). 

It should be noted that the Ising model with boundary magnetic field has been studied 
much, especially on the lattice and on the half-line. We refer the reader to |I2|, [23l [26l [27] 
and the references in them for further information. 

To recognize conformal symmetry and to identify the representation content it is suf- 
ficient to look at the ratios of the energy gaps; therefore one often considers normalized 
spectra which are obtained by subtracting the ground state energy and dividing by the low- 
est energy gap. The normalized exact and TCSA spectra for the fiows ((31), ([4]) as a function 
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of the logarithm of h can be seen in Figure [H An interesting feature of these TCSA spec- 
tra is that they appear to correspond to the flows (1/2, 1/16) (1/2, 1/2) — > (1/2, 1/16) 
and (1/2,1/16) (1/2,0) (1/2,1/16), i.e. second flows appear to be present after 
the normal flows. Flows in other minimal models also show this behaviour [201 ES]. One 
application of the picture ^ could be the explanation of this phenomenon. 

Following Watts' proposal based on the look of the TCSA spectra shown in Figure [8] 
we assume that only the flrst two terms are nonzero in ([T]): 

H'^h^ric) = So{h,nc)Ho + Si{h,nc)Hi . (5) 

In summary, we look for answers for the following questions: 1. Does the spectrum of 
(l5|) agree with the TCSA spectrum in some approximation with a suitable choice of the 
functions Sq and Si? 2. How can we explain the 'second' flows in the TCSA spectra? 

We determine the exact spectrum of (l2|) using an essentially known (see e.g. [23l [26l [271 
EH EH [30]) quantum held theoretic representation of this operator. In this representation 
the operator ([21) is a quadratic expression of fermionic flelds. We extract the spectrum 
from the (quantum) field equations, which are linear. 

Concerning the TCSA Hamiltonian operator, the main difficulty we encounter is that 
it is hard to handle the TCSA spectra analytically even if the non-truncated model is 
exactly solvable. Therefore, hoping that we can gain some insight by looking at a similar 
but exactly solvable truncation method, we tried another truncation method which we 
call mode truncation. The mode truncated model can be solved exactly, but it turns out, 
rather unexpectedly, that the behaviour of the spectrum for large values of h is different 
from the behaviour of the TCSA spectrum, namely the qualitative behaviour of the mode 
truncated spectrum is very similar to that of the exact spectrum; the second flows are not 
present. This, besides leaving the second question open, raises the problem of flnding the 
possible behaviours for large values of h and their dependence on the truncation method, 
and whether the mode truncation method can be generalized to other models. 

We also apply the Rayleigh-Schrodinger perturbation theory to verify the validity of 
([5]) for both the TCSA and mode truncation methods. 

A further calculation that we do is a numerical comparison of the exact and TCSA 
spectra. We perform the same calculation for the mode truncation scheme as well. 

Our secondary aim in this paper, which is partially independent of the problem of the 
TCSA approximation, is to provide a detailed description of the quantum held theoretic 
representation of Although this was studied in the literature (e.g. in [28] and especially 
in [26] (at flnite temperature), and on the half-line in several papers), we give a description 
in which the perturbed Hamiltonian operator structure, on which the TCSA is based, is 
emphasized and observed consistently. In a calculation using the Bethe-Yang equations 
(which we also present in the appendix) or in the approach of |26| where boundary condi- 
tions play central role, the link with the TCSA formulation would not be entirely obvious. 
Furthermore, the formulation we present is also suitable for Rayleigh-Schrodinger pertur- 
bation theory and for the treatment of the mode truncated version. We do not consider the 



5 



classical level of the field theoretic model, mainly because it is irrelevant to our problem. 
We remark that the same model with massive unperturbed part (see e.g. [12]) could be 
studied along the same lines. The field theoretic approach also raises the problem of defin- 
ing distributions (or similar objects) on a closed interval. We do not know of a systematic 
exposition of this subject (neither for a closed interval nor for the half-line), although it 
would be needed for boundary field theory. In this paper we use distributions on closed 
intervals, restricting to the most necessary formulae only. These formulae are collected in 
Appendix [Al 

We remark that our TCSA program relies on the conformal transformation properties 
only and does not make use of the representation mentioned above. 
The contents of the paper are the following: 

Section [2] contains the field theoretical description of the model (l2|), in particular the 
calculation of the exact spectrum. Further details concerning the boundary conditions, 
the normalization of interacting creation and annihilation operators, the relation between 
free and interacting creation and annihilation operators, matrix elements of the interacting 
fields, identities for the Dirac delta on the strip and the expression of eigenstates in terms 
of the unperturbed eigenstates are deferred to Appendix [Bl 

In Section [3] we propose a description of (l2|) as a perturbation of the h ±oo limiting 
case, which cannot be found in the literature. An interesting feature of this description, 
which we call reverse description, is that the perturbing operator is non-relevant. We 
calculate the exact spectrum in a similar way as in the case of the standard description 
mentioned above. This reverse description is motivated by the presence of the second 
(reverse) fiows in the TCSA spectra. 

In Section [4] we describe the mode truncated model and calculate its spectrum. This is 
done along the same lines as in the non-truncated case. 

Section [5] contains the perturbative results for the sq and si functions. 

Section [6] contains the results of the numerical test of the approximation by (JH]) for the 
TCSA and the mode truncation schemes. 

The paper contains further appendices beyond those mentioned above. 

In Appendix [Cl we describe the spectrum using the Bethe-Yang equations. We find that 
they give the exact result in this case. 

Appendix [P] contains the power series for the exact, TCSA and mode truncated energy 
levels up to third order in h which we obtained by means of the Rayleigh-Schrodinger 
perturbation theory. 

Conclusions and some directions for future work are given in Section [71 

2 The exact spectrum 

In this section we present a quantum field theoretic description of the model ([2]) and 
calculate its energy spectrum. 
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2.1 The free model 

The defining constituents of the unperturbed model are the following: two fermion fields 
t) and ^'2(0;,^) and a fermionic operator A2{t) with the anticommutators 

{^i(a;, t), -^^iy, t)} = {^2(2;, t), ^^2(1/, t)} = AL6{x - y) (6) 

t), ii/2(l/, t)} = -AL{5{x + y) + 5{x + y- 2L)] (7) 

{A2(t),^i(x,t)} = {A2(i),*2(x,t)} = (8) 

{^(i),^(i)} = 2 (9) 

and the relations 

t)^ = t) *2(X, t)^ = *2(X, ^2(t)^ = ^2(i) , (10) 

the Hamiltonian operator 

Ho^-^l 0)^-*i(^> ^^^izl 0) , (11) 

the equations of motion 

^^A,{t) = [iHo,A,{t)] = (12) 
dt-^i{x,t) = [iHo,-^i{x,t)] 

= -d,^i{x, t) + ^m-x) + e{x - L)][-d,^2ix, t) + d,^i{x, t)] 

+ ^[-^1(0, t) - ^2(0,t)]5(a;) + l[^i{L,t) + ML,tMx - L) , (13) 
dt'^2{x,t) = [ii/o,*2(2;,i)] 

= a,*2(x,i) + ^M-x) + e(x - L)][-d,Mx,t) + d,^i{x,t)] 

+ ^[*i(0, t) + *2(0, t)]d{x) + ^[-*i(L, t) - *2(i^, t)]5(a: - L) . (14) 

The fermion fields, which are one-component real fermion fields with zero mass, have the 
following expansion: 

^i{x,t)= J2 \/2[a(A;)e^^(*-^) + at(A;)e-''=(*-^)] + A (15) 

fcgf Z, fc>0 

^2(x,t)= V2[-a(A;)e^*^(*+^) - a1-(fc)e-^'=(*+")] - (16) 

fcef Z, A;>0 
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a{k)^ = a{-k) A\ = (17) 

{aiki), a{k2)} = 5k„^k, {aik), A^} = (18) 

{A2,A,} = {A,, A,} = 2 (19) 

{A2,a{k)} = 0. (20) 

The fermion fields and A2 are dimensionless. 

A unitary representation for the above operator algebra is defined by the following 
formulae: an orthonormal basis of the Hilbert space is 

a{ki)a{k2)a{k3) . . . a{kn)\u) a{ki)a{k2)a{ks) . . . a{kn)\v) , (21) 

where ki > k2 > ks > ■ ■ ■ > kn > 0, ki E ttZ/L and n > 0, 

Aia{ki)a{k2)a{k3) . . . a{kn)\u) = (-l)"a(A;i)a(A;2)a(A;3) • • • a{kn)\v) (22) 

Aia{ki)a{k2)a{ks) . . . a(A;„)|f ) = (-l)"a(A;i)a(A;2)a(/c3) • • • ci{K)\u) (23) 

^2a(A;i)a(A;2)a(/c3) . . . a{kn)\u) = (-l)"(-i)a(fci)a(fc2)a(/^3) • • • a{kn)\v) (24) 

A2a{ki)a{k2)a{k3) . . . a{kn)\v) = {-iy\+i)a{ki)a{k2)a{k3) . . . a{kn)\u) . (25) 



The condition ki > k2 > k^ > ■ ■ ■ > kn > 0, ki E vrZ/L and n > also applies below if 
not stated otherwise. a{k)\u) = 0, a{k)\v) = if A; < 0. The energy eigenvalue of a basis 
vector a{ki)a{k2)a{k3) . . . a{kn)\u) is ^^=1 ki, and the same applies if m is replaced by v. 
The following boundary conditions are satisfied by the fermion fields: 

(Ei|^i(0,t) + ^2(0,t)|^2) = (Ei|^i(L,t) + <il2mt)\E2) = , (26) 

where \Ei) and {£2) are arbitrary energy eigenstates. 

A representation of the Virasoro algebra can be defined on the Hilbert space in the 
following way (which is well known essentially in conformal field theory; see [3TI 132] ): 

LN = ^y] -ka{-k)a{k - ^) (27) 

^0 = ^ E ^ <-k)a{k) :] + ^ , (28) 

where N = ±1,±2,..., a(0) = (l/-\/2)Ai, :: denotes the normal ordering for fermionic 
creation and annihilation operators, : a{—k)a{k) := a{—k)a{k) if < 0, : a{—k)a{k) := 
—a{k)a{—k) if /c > 0. It can be verified that in the above representation of the a{k) 
and A2 operators L^, N E Z satisfy the usual relations of the generators of the Virasoro 
algebra with central charge c = 1/2. A2 commutes with the L^- \v) and \u) are highest 
weight states with weight 1/16 and the Hilbert space decomposes into two copies of the 
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M(c = 1/2, h = 1/16) unitary highest weight representation of the Virasoro algebra. The 
invariant subspace belonging to u is spanned by the vectors a{ki)a{k2) ■ ■ ■ a{kn)\u) with 
n even and a{ki)a{k2) ■ ■ .a{kn)\v) with n odd. The invariant subspace belonging to v is 
spanned by the vectors a{ki)a{k2) . . . a{kn)\u) with n odd and a{ki)a{k:2) ■ ■ ■ a{kn)\v) with 
n even. These subspaces will be called u and v sectors. The relation between Hq and Lq is 

= - 1^ . (29) 

The fields *i and *2 can be written as ^1(2;) = V^J2n& o.{n)z'', *2(^) = V^J2n& -aC"-)^", 
where a(n) = a{k) , k = mr / L, z = exp[ij(t — x)], z = exp[i^(t + ,x)]. For e{z) = '^l{z)/^/z 
and e{z) = ^2(^)/Vi we have [Ljv,e(z)] = ^^+^^(^) + i(iV + l)z^e(z) and [LN,e{z)] = 
z^~^^^^{z) + 7^{N + l)^^e(^), i.e. e{z) and e{z) are chiral Virasoro primary fields of weight 
1/2. The same relations apply to the fields A2e and A2e, in particular A2e and A2e are also 
chiral primary fields of weight 1/2. A2e and ^426 have zero matrix elements between the u 
and V sector, whereas e{z) and e{z) have zero matrix elements within the u and v sectors. 
In the above equations the domains of z and z are extended to the whole complex plane. 

The operator A2 is an auxiliary operator and if it is omitted, then it is possible to 
represent the fields so that the Hilbert space is a single Virasoro module (1/2, 1/16). It is 
the next section where the presence of A2 will be really useful. 

2.2 The perturbed model 

The Hamiltonian operator is 

H^—^J 0)d,^i{x, 0)dx +-^J ^2{x, 0)d,^2{x, 0) dx 

+ hiA2{0)[^2{L, 0) - *i(L, 0)] , (30) 

where /i is a coupling constant of dimension mass. The perturbing term 

hHi = hiA2{0)[^2{L, 0) - *i(L, 0)] (31) 

has zero matrix elements between vectors belonging to different sectors, which means that 
H can be restricted to the u and v sectors separately. These restrictions are denoted 
by H\u and Hj is also a primary boundary field of weight 1/2 with respect to 

the Virasoro algebra representation defined in the previous section taken at t = 0. The 
matrix elements of Hj, i.e. of Hi\u and Hi\y, are uniquely determined by this property 
and by the values of the matrix elements {u\Hi\u) and {v\Hi\v). It is easy to verify that 
2 = {u\Hi\u) — —{v\Hi\v), and there exists an intertwiner Y of the Virasoro algebra 
representations on the u and v sectors so that Yu = v, YHo\uY~^ = HqI^. This also 
implies that YHi\uY~^ = —Hi\y. This means finally that we can restrict to < /i < 00, 
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and the u sector and H\u will correspond to the h > case of ^ and to ((HI), the v 
sector and will correspond to the h < case of ([2]) and to dH). For < h < oo the 
operator ( l30ll describes in the two sectors the two flows mentioned in the Introduction. The 
precise relation between the /i in ([2]) in the Introduction and the h in f l30ll is the following: 
hintrod = 2L^/'^h in the u sector, hintrod = — 2L^/^/i in the v sector. Further on we shall 
assume that < h < oo. 

The eigenvalues of the Hamiltonian operator (l30l) are ultraviolet divergent in perturba- 
tion theory; the divergence can be removed by adding a term ch^I with appropriate value 
of the logarithmically divergent (as a function of the cutoff energy) coefficient c. This 
means that the differences of the eigenvalues of H are not ultraviolet divergent. We shall 
assume that the ground state energy is set to zero, and the c/i^J term will not be written 
explicitly. 

The equations of motion are obtained by adding the following terms to the right-hand 
side of the unperturbed equations of motion (fT2i l. (fT3ll . ( fT4l ): 

[ihHi, t)] = 8LhA2{t)6{x - L) (32) 
[ihHi, ^2(2;, t)] = -8LhA2it)S{x - L) (33) 
[ihH:, A2{t)] = 2h{m2{L, t) - ^i(L, t)) . (34) 

The initial condition for the fields \E'i(x, t), \E'i(x,)f:) and A2(t) is that at t = they are 
equal to the unperturbed fields (of the previous section). The above terms are linear in 
the fermion fields and A2, so the equations of motion for the perturbed theory are linear 
and by sandwiching these equations between energy eigenstates we get a system of three 
first-order differential equations for the expectation values of the fields and A2. These 



expectation values can be assumed to take the form 

(Ei|^i(x,t)|E2) = -e^^(*-^) - Q{x - L)Ci{k)e'^' (35) 

(El 1^2(3^,^)1^2) = e''=(*+^) - e{x - L)C2{ky^' (36) 

{Er\A2{t)\E2) = C,{ky^\ (37) 



where Ci(A;), C2(A;), C-i{k) are finite constants, \Ei) and |i?2) are eigenstates of H and 
k = El — E2. \Ei) and I-E2) are not necessarily normalized to 1 here. Substituting ([35l) - 
f l37ll into the equations of motion we get algebraic equations for Ci{k), C2{k), Cs{k), which 



have the following solution: 

kLtaii{kL) = WL'^h^ (38) 

^i(A;)(x,t) = {Ei\'^i{x,t)\E2) = -e^^(*'^) - 6(0; - L)i sin(fcL)e^^* (39) 

Mk){x,t) = {Ei\^2ix,t)\E2) = e''=(*+") - e(x -L)isin(A;L)e''=* (40) 

a2m) = {EM2it)\E2) =_i^i^e^H_ ^4^^ 
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Equation (13811 is the formula that determines the possible values of at a given value of h 
and L and so the spectrum of H up to an undetermined additive overall constant. 

The assumption (l35l ). ( l36l ) is motivated by the fact that the equations of motion for 
X G (0,L) are {dt + dx)'^i{x,t) = 0, {dt — dx)'^2ix,t) = 0. We remark that the equations of 
motion could also be solved directly without making any assumptions on the form of the 
expectation values. 

Introducing the notation 

n{k) = {Mk),Mk),a2ik)) (42) 
the mode expansion of \l/2, ^2) is 

{^,,^2,A2) = Y,bikMk) , (43) 

kes 

where the b{k) are creation/annihilation operators and the summation is done over the set 
S of all real solutions of (l38l) . The b{k) satisfy the relations (see Appendix [B] for details) 

WkM>=M=K... ,,,^'%,,^^ (44) 

and b{ky = b{-k). 

The Hilbert space is spanned by the orthogonal eigenstates b{ki)b{k2)b{k^) . . . 6(/c„)|0/i) 
where fci > A;2 > ■ ■ ■ > > 0, |0/i) is the ground state, which is unique, and b{k)\Qh) = if 
< 0. The eigenvalues of these states are ^^=1 ki- The eigenvector 6(^1)6(^2) • • • b{kn)\Qh) 
belongs to the v sector if n is even and to the u sector if n is odd. The first few energy 
gaps (i.e. energies relative to the lowest energy) within the two sectors are shown in Figure 
[7] as functions of ln(/i) with L = 1. 

In the /i : — > 00 limit 

k{n, h) ^ fc(n, 0) + -- = — + -- , (45) 

where k{n,h) is the n-th nonnegative root of ([38l) as a function of k{n,0) = mr/L, 
n = 0,1,2, ... ; and {b{ki), b{kj)} — » 26ki+kj,o- It can be verified that in the /i — > 00 limit 
the (c = 1/2, /i = 0) representation of the Virasoro algebra can be introduced in the v sector 
and the (c = 1/2, /i = 1/2) representation can be introduced in the u sector. One can write 
expressions (which are well known essentially, see [Slilffij) for the generators in terms of the 
b{k) similar to ((27]), ((28]). Therefore h : ^ 00 corresponds to the (1/2, 1/16) (1/2,0) 
flow in the v sector and to the (1/2, 1/16) (1/2, 1/2) flow in the u sector. We remark 
that the easiest way to determine which representations are realized in the h ^ 00 limit is 
by counting the degeneracies of the first few energy levels (separately in the two sectors). 
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3 Reverse description 



In this section we propose the description of the model (l30ll as a perturbation of its /i — > oo 
limit. This is motivated by the presence of the second (reverse) flows in the TCSA spectra 
(Figure [HI). It should be noted that the meaning of the notation Hq or Hj etc differs from 
that in Section O The precise correspondence between the quantities in this section and in 
the previous section will be given explicitly for the coupling constant and for the spectrum. 

3.1 The free model 

The fundamental objects of the model at /i = oo are two one-component real massless 
fermion fields $i(x,t), $2(2;, t) with the anticommutation relations 



and reality property $i(a;, t)^ = t), $2(2;, t)^ = $2(2;, t). $1 and $2 are dimensionless. 
The Hamiltonian operator is 



{<l>i(x, t), <l>i(y, t)} = {$2(x, t), <l>2(y, t)} = 4L6{x - y) 
{$i(a;, t), $2(1/, t)} = -AL[5{x + y) - 5{x + y - 2L)] 



(46) 
(47) 




dx <l'i(x,0)a^$i(a;,0) + 




L 



dx $2(0;, 0)9^$2(a;,0) . 



(48) 



The equations of motion are 



dt^i{x,t) = [iHo,<l>i{x,t)] = -d,<!>i{x,t) 

+^Six - L)[$i(L,t) - $2(i^,t)] + ^5(x)[-$i(0,t) - $2(0, t)] (49) 

+^0(-x)[9,$i(x,t) - 6',$2(x,t)] + ^0(x - L)[5,$i(x,t) + d,<!>2{x,t)] 



+ ^S{x - L)[$i(L,t) - ML,t)] + ^6{x)[<^^{0,t) + $2(0, t)] (50) 
+ le(-x)[9,$i(x,t)-9,$2(x,t)] + ^0(x-L)[-9,.$i(x,t)-9,$2(x,t)] . 



The fermion fields have the following mode expansion: 





(53) 
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{a{ki),a\k2)} = Sk^^k2 a{ky = a{-k) . (54) 

An orthonormal basis for the Hilbert space is formed by the vectors a{ki)a{k2) ■ ■ ■ a(A;„)|0), 
where h > 0, k e ttZ/L + 7r/(2L), n > 0. a{k)\0) = if < 0. The eigenvalue of 
the eigenvector a{ki)a{k2) ■ ■ ■ a{kn)\0) is ^"^^ ki. The fermion fields satisfy the following 
boundary conditions: 

(Ei|<l>i(0,t) + $2(0,t)|E2) = (^i|$i(L,t)-<l>2(L,t)|E2) = 0, (55) 

where \Ei), \E2) are eigenstates of i^o- 

One can define a representation of the Virasoro algebra on the Hilbert space in the same 
way as in Section O (see also [3ll[32]). This representation is (1/2,0) © (1/2,1/2). The 
fields $i(x, t) and $2(2;, t) can be converted to weight 1/2 primary fields by multiplying 
them by a suitable simple exponential factor. 

3.2 The perturbed model 

The perturbed Hamilton operator is 

H = --^ [ dx <l>i(x, 0)9,$i(x, 0) + -V /" dx $2(x, 0)9,.$2(x, 0) 
oJ^ Jo oL Jq 

+ (7[i($i + $2)(L,0)lim9,.($2-$i)(x,0)] , (56) 

where g is a dimensionless coupling constant. In the same way as in Section 12.21 the 
perturbing term has zero matrix elements between vectors belonging to different irreducible 
representations. 

The limit prescription in the perturbing term is important, and the limit lim^.^^ should 
be taken at the end of any calculation. It should be assumed that limx^L6{x — L) = 0, 
for example, and similarly for the derivatives of 6{x — L). 

It should be noted that the only primary field in the or in the 1/2 representation is 
the identity operator, all other fields are descendant and non-relevant fields. 

The equations of motion are obtained by adding the following terms to the right-hand 
side of the unperturbed equations of motion f l49l ). (l50ll : 

[igHi, $i(x,t)] = g8L5{x - L) lim 9^.($2 - $i)(x,t) (57) 
[igHj,<!>2{x,t)]= g8L6ix-L) lim 9,($2 - $i)(x, t) (58) 

X—rL 

where 

i/,= [i(<l>i + <l>2)(L,0)lim9,(<l>2-$i)(a:,0)] . (59) 
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Similar steps to those in Section [2?2] can now be taken to obtain t)|i?2) and 

{Ei\^2{x,t)\E2) . In the same way as in Section [2?2l the forms 

(Ei|$i(s,t)|E2) =e''=(*-") + 0(x-L)Z}i(fc)e''=* (60) 
{Ei\<^2{x,t)\E2) = -e^'=(*+^) + Q{x - L)D2{k)e"'' (61) 

can be assumed, where Di{k) and D2{k) are finite constants, \Ei) and l-E'2) are eigenstates 
ofH and k = Ei - E2. 

Solving the equations of motion for Di{k), D2{k) we get D2{k) = —Di{k) = cos{kL) 
and 

(kL) tan(A;L) = — . (62) 
16g 

Equation (l62D is the formula that determines the spectrum of H up to an overall additive 
constant. The eigenvalues of H are ki + k2 + ■ ■ ■ + kn, where n > 0, ki > 0, ki ^ kj if 
i ^ j, the ki are real roots of (l62l) and the lowest eigenvalue is assumed to be set to zero 
by adding a constant which is not written explicitly. The substitution 

' = T^hh- ^^^^ 

converts (l62ll into (l38l l. Thus g : Q ^ —00 corresponds to the flow (0 1/16) © (1/2 
1/16). 

In perturbation theory there are divergences if we take x = L in Hj at the beginning. 
However, if one allows x to take general values, then in Rayleigh-Schrodinger perturbation 
theory for the differences of the energy levels one can expect to get sums at any fixed order 
which are possible to evaluate. We expect that the evaluation yields, besides non-singular 
parts, step functions, 6{x — L) and its derivatives, and the result remains finite after the 
a; — > L limit. 



4 Exact spectrum in the Mode Truncated version 
4.1 The free model 

Let called the truncation level, be a positive integer. The mode truncated version of 
the free model described in Section 12.11 is the following: 

{^i(x,t),^i(y,t)} = 2[l + 2 cos(A;(x - y))] (64) 

fcef {l...n4 

{^2(x,t),^2(2/,t)} = 2[l + 2 cos(A;(x - y))] (65) 

fcef {l...n4 

{$i(x, t), $2(1/, t)} = -2[1 + 2 Y + y))\ ' (66) 

fcef {l...nc} 
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feef {l...nc} 

^2(0;, t)= V2[-a(A;)e^^(*+") - a+(fc)e-'^(*+")] - A, . 

fcef {l...n4 

Equations (fT7l)-(l20ll apply unchanged. The Hamiltonian operator is 

Ho= k[a{k)a'^{k)] . 

fcef{l...n,} 

The fields satisfy the following equations of motion: 

^^A2{t) = [iHo,A2{t)] =0. 
The Hilbert space and the energy eigenstates are similar to those in Section [2 

k e -j^i-^c, -Uc + 1, . . . , nc - 1, nj 
applies instead of A; G TxTLjL. The Hilbert space is 2 x T^" dimensional. 

4.2 The perturbed model 

The perturbed Hamiltonian operator is i^o + hHi-, where 

i/, = iA2(0)[^2(A0)-^i(L,0)] . 
The equations of motion are 

di^^{x,t) = [i{Ho + hHj),^,{x,t)] = -d,^^{x,t)-hA2{t)C,{x) 

dt<il2ix,t) = [i{Ho + hHj), ^2ix,t)] = d,^2ix,t) - hA2{t)C2{x) 

^^A2{t) = \x[E^^hEi),A2{t)\ =2/i(^2(At)-^i(At)) , 

where 

= -4[l + 2 Y cos(A;(x + L))] 

fcef{l...ne} 

C2{x) = {^2{x,t),<il2{L,t) - ^i(L,t)} = -Ci{x) . 
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These equations are linear as in the non-truncated case, so sandwiching them between 
energy eigenstates gives a system of three first-order linear partial diff'erential equations 
for the expectation values. The analogue of (l38l l can be obtained from these equations in 
the following way: we can eliminate A2: 

d^^,{x,t) = -d,t^i{x,t) - 2h\^2{L,t) - ^i{L,t))Ci{x) (80) 
d^^2ix,t) = 9,i*i(x,t) - 2h\^2mt) - ^i(L,t))C2(x) (81) 

and introduce the functions /i(x), f2{x): 

(Ei|^i(x,t)|E2) = /i(x)e^^* {E,\^2{x,tm) = f2{x)e'^' , (82) 

where k = Ei — E2 and the dependence of /i and /2 on k is not denoted explicitly. 
Sandwiching ( l80l) and ( ISTll between \Ei) and \E2) gives two first-order inhomogeneous 
linear ordinary differential equations (with the deviation that the inhomogeneity depends 
on the unknown functions). We also have the boundary conditions /i(0) = — /2(0), fi{L) = 
—f2{L). Elementary calculation yields that a solution of the differential equations and the 
boundary conditions exists if and only if k satisfies the equation 



fc2 + 5Z 



1.2 Z-^ ^2 _ 1.2 

fcoef {l...nc} " 



(83) 



which is the analogue of (l38l) and determines the energy of the modes as functions of h. 

fl83l) as an algebraic equation for k has finitely many real roots, and if is a root, then 
—k is a root as well. All real roots converge to finite values as /t — *^ 00 except for the pair 
with the largest absolute value. This pair of roots diverges linearly as — cx). This pair 
has the largest absolute value already at h = 0. A consequence of this behaviour is that 
the lower half of the spectrum, namely those states which do not contain the mode with 
the highest energy, remains finite as /i — 00, whereas the higher half of the spectrum, i.e. 
the states which contain the mode with the highest energy, diverges linearly as /i — > 00 
with a common slope. Here it is assumed that the ground state energy is set to zero. In 
the subsequent sections we shall consider the lower half of the spectrum. The look of this 
half as a function of \ia{h) is very similar to that of the exact spectrum shown in Figure [71 

Applying the formula 

^ 2k^ _ k/2 
+ 2^k2_ 4^2^2 - tan(A;/2) ^ 



n=l 



we can easily verify that the limit of (l83ll as ric ^ 00 is ([38 
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5 Perturbative results 



The functions So{h,nc) and Si{h,nc) in = SQ(h,nc)Ho + Si{h,nc)Hj are determined by 
the condition that the differences of those eigenvalues of that are low compared to the 
truncation level should be equal to those of the truncated Hamiltonian operator H^ijic). 
H\nc) is, in particular, the TCSA Hamiltonian operator H'^'-^^^^ric), or the Hamiltonian 
operator H^'^'^ijic) of the mode truncated model. This condition, which we shall call 
renormalization condition, applies separately and independently within the u and v sectors, 
and we have in fact a pair Sg, for the u sector and another pair Sq, s\ for the v sector. 

The renormalization condition is a very strong condition on Sq and Si and generally we 
cannot expect that it can be satisfied. It is possible, however, that it can be satisfied in 
certain approximations. 



5.1 Mode Truncation scheme 

Using (11561) in Appendix [D] we can obtain the following results: 

The renormalization conditions have a solution if the eigenfunctions are expanded into 
a power series in h and terms that are higher order than 3 are omitted: 

so{h,n^) = l + 0{h^) si{h,n^) = h + ]^{S-S{n^))L^h^ + 0{h^) , (85) 



where 



°° 32 . 32 



n=l n=l 



This solution applies to both the u and v sectors and it is exact in ric. We remark that 

32 



S-S{n,) = ^ + 0{l/nl) . (87) 



In the MT scheme we can obtain another result that is non-perturbative in h: doing 
power series expansion in l/ric we obtain the formula 

Omitting the terms which are second or higher order in 1/nc, equation (l83l) takes the form 
of ([38]) if 

h 1 

Ks= , ^= = h + lQh^ — — + 0{l/nl) (89) 

is introduced: 

kLi&n{kL) = lQL^hl^ . (90) 
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This means that rescaling by 

so{h,n,) = l + 0{l/nl) si{h,n,) = h + IQh^ — — + 0{l/nl) (91) 

improves the convergence of the mode truncated energy gaps to the exact energy gaps 
from order l/ric to order 1/n^ (at least), i.e. the difference between the energy gaps of 
so(h,nc)HQ + si{h,nc)Hi and the energy gaps of {Hq + hHj)^'^ tends to zero as 1/ra^ (at 
least) for any fixed finite value of h, whereas the difference between the energy gaps of 
Hq + hHj and {Hq + hHiY'^"^ tends to zero as l/ric- 



5.2 TCS scheme 

Using 0156p in Appendix [D] we obtain the following results: 

The renormalization conditions have a solution if the eigenfunctions are expanded into 
a power series in h and in l/rZc and terms that are of order higher than 3 in /i and 1 in 
l/ric are omitted: 

SQ{h,n^) = l+xih'^ + yih^ + 0{h^) (92) 
si{h,n^) = h + y2h'^ + X2h^ + 0{h^) (93) 

xi = ^ + 0(lK) yr = Q + 0{l/nl) (94) 

^2 = \{S- S{n,))L^ + 0{l/nl) 1/2 = + 0{l/nl) . (95) 

This solution applies to both sectors. If terms of order l/n1 are also taken into consider- 
ation, then the renormalization conditions do not have a solution. A significant difference 
between the TCS and mode truncation schemes is that the value of the coefficient xi is 
zero in the mode truncation scheme but non-zero in the TCS scheme. 



6 Numerical results 

In the numerical calculations described in this section the value of L is set equal to 1. 
This does not affect the generality of the results. In the calculations we used the same 
normalizations as in Section [2l 

so{h,nc) and si{h,nc) are calculated from the lowest three energy levels and are given 
by the following formulae: 
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where the superscriptless quantities are the non-truncated ones and those with the * super- 
script are obtained by a truncation. Eq, Ei, E2 denote the three lowest energy eigenvalues. 



6.1 Mode Truncation Scheme 

Figure [T] shows the exact and mode truncated spectra as a function of the logarithm of the 
coupling constant. The truncation level is nc = 9, and the dimension of the Hilbert space 
is 512 in each sector. It is remarkable that there is a good qualitative agreement between 
the mode truncated and exact spectra for all values of h. 

Figure [2] shows the same spectra, but the lowest gap is normalized to one, i.e. the 
functions are shown. It is remarkable that the agreement between the exact 

and mode truncated spectra looks considerably better than in the case of not normalized 
spectra. 

Figures [3H5] show the functions so{h), Si{h), si{h) / so{h) determined by the lowest three 
energy levels in the v sector via the formulae (i96l ). ( l97l l in various ranges. Figures [4] and 
[5] also show the curves given by (l85ll on the left-hand side (red/grey line). So{h) remains 
close to 1 and for large values of h it tends to a constant which can be expected to converge 
to 1 as nc 00. si{h) also tends to a constant for large values of h which can be expected 
to increase to infinity as nc ^ 00. The behaviours of Si{h) / so{h) and si{h) are similar. 

Figure Oa shows the normalized mode truncated spectrum and the normalized exact 
spectrum rescaled by so{h) and si{h) (i.e. the normalized spectrum of H^) in the v sector. 
No difference between the two is visible. In Table [T] values of the fifth normalized energy 
gap fc(i'fe)^fc|o'fe) of the v sector are listed: the values in the non-truncated case are listed 
in the first column, the values in the mode truncated case are listed in the second column 
and the rescaled non-truncated values are listed in the third column (which would be the 
same as the values in the second column if the renormalization condition could be satisfied 
exactly) . 

Figure [6lb shows the normalized mode truncated spectrum and the normalized exact 
spectrum rescaled by sq^H) and si{h) in the u sector. The so(/i), si{h) obtained in the v 
sector were used for the rescaling, which corresponds to the assumption that so^h) and si{h) 
are the same for both sectors. The difference between the mode truncated and rescaled 
exact spectra is not visible in the figure. 

We also see from Table[2]in which values of the fourth normalized energy gap fc(i'fe)Ifc(o'fe) 
of the u sector are listed that the rescaling together with the above assumption works well. 

We have not tried to calculate sq and si for the u sector because (h) is not 

invertible in this case. One way to circumvent this difficulty would be to use other energy 
levels Ei,Ej,Ek for which ^^^(h) is invertible. 
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Figure 2: Exact (dashed lines) and mode truncated (solid lines) normalized spectra in 
V and u sectors respectively as a function of ln(/i) at truncation level ric — 9 
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Figure 3: The function SQ{h) for the v sector in the ranges /i G [0, 3], G [0, 20], G [0, 400] 
and Sq G [0.95, 1.05] at truncation level = 9 
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Figure 4: The function Si{h) for the v sector in the ranges he [0, 2], G [0, 20], /i G [0, 400] 
and Si G [0, 2] at truncation level ric = 9 
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Figure 5: The function Si{h)/so{h) for the v sector in the ranges h G [0,2], h G [0,20], 
h G [0,400] and Si/sq G [0,2] at truncation level ric = 9 
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Table 1: The normalized energy gap fc(i'/tj_|_fc(o'ftj in the v sector: exact, MT {ric = 9) and 
rescaled exact values. 
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Table 2: The normalized energy gap k{ifJj_}}^Qh) the u sector: exact, MT {ric — 9) and 
rescaled exact values. 



ln(/i) Exact MT Rescaled exact 
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6.2 TCS scheme 



Figure [7] shows the exact and TCSA spectra as a function of the logarithm of the coupling 
constant. The truncation level is ric = 14, and the dimension of the Hilbert space is 110 in 
each sector. It is remarkable that there is strong deviation between the TCSA and exact 
spectra for large values of h. The behaviour of the TCSA energy gaps is Ei{h) — Eq^K) oc h 
for large values of h. 

Figure [8] shows the same spectra, but the lowest gap is normalized to 1, i.e. the functions 
cire shown. It is remarkable that the agreement between the exact and TCSA 

spectra looks better than in the case of not normalized spectra. The functions Et(h)-^Eo{h) 
have finite limit as ^ oo and the degeneracy pattern in this limit appears to correspond 
to the c = 1/2, h = 1/16 representation of the Virasoro algebra. The correspondence 
improves as ric is increased (this improvement is not illustrated). At any fixed finite value 
of h, however, the TCSA data are expected to converge to the exact values as nc oo. 

Figures iQlfTTl show the functions so{h), si{h), si{h)/so{h) in various ranges calculated 
in the same way as in the mode truncated case. The figures also show the curves given 
by (l92l) . (l93ll on the left-hand side (red/grey line). It is remarkable that So(/i) oc h for 
large values of h. Calculations at other values of ric show that the slope of So{h) decreases 
as ric is increased and it can be expected to converge to as ric —> cxd. si{h) appears to 
tend to a constant for moderately large values of h. Calculations at other values of ric 
show that this constant increases as ric is increased and it can be expected to converge to 
infinity as oo. For large values of h, Si{h) decreases. Si{h) / so{h) reaches a maximum 
at /i ~ 1.6 and then decreases to zero. Calculations at other values of ric show that the 
maximum value and the value of h where it is reached increase as ric is increased and it 
can be expected that both values converge to infinity as nc oo. 

Figure [T2I a shows the normalized TCSA spectrum and the normalized exact spectrum 
rescaled by So{h) and Si{h) (i.e. the normalized spectrum of H^). They show good quali- 
tative agreement. 

Values of the fifth normalized energy gap fc(i'fe)^fc(o'fe) of the v sector are listed in Table 
[3] as in the mode truncated case. These data show that the rescaling significantly improves 
the agreement between the TCSA and exact spectra (which is especially noticeable if 
ln(/t) > -2). 

Figure [T2lb shows the normalized TCSA spectrum and the normalized exact spectrum 
rescaled by so{h) and si{h) in the u sector. The so{h), si{h) functions obtained in the v 
sector were used for the rescaling as in the mode truncated case. Table [4] shows values of 
the fourth normalized energy gap of the u sector. We see from the table that 

the assumption that so{h) and si{h) are the same in both sectors does not give very good 
result in this case, although the situation might become better at higher truncation levels. 
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Figure 7: Exact (dashed lines) and TCSA (solid lines) energy gaps {Ei — Eq) in the v and 
u sectors respectively as a function of ln(/i) at truncation level Uc = 14 




Figure 8: Exact (dashed lines) and TCSA (solid lines) normalized spectra in the v and u 
sectors respectively as a function of ln(/i) at truncation level ric — 14 
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Figure 10: The function Si{h) for the v sector in the ranges h G [0,1.75], Si G [0,1.7]; 
h G [0,20], si G [0,1.5]; h G [0,100], Si G [0,1.5]; h G [0,400], Si G [0,1.5] at truncation 
level rir = 14 
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Figure 11: The function Si{h)/so{h) for the v sector in the ranges h G [0,1.75], Si/sq G 
[0, 1.75]; h e [0, 10], Si/sq e [0,1]; he [0,400], Si/sq G [0, 1] at truncation level = 14 
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Table 3: The normalized energy gap fc(i'fej_|_fc(o'fcj in the v sector: exact, TCSA ( 
and rescaled exact values. 

\n{h) Exact TCSA Rescaled exact 
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Table 4: The normalized energy gap fc(i'/tj_fc(o'/tj in the u sector: exact, TCSA ( 
and rescaled exact values. 
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3. 


,046206 


3. 


,046198 


-3 


3. 
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3. 
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7 Discussion 



We have investigated the validity of the approach ([5]) for the description of truncation effects 
in TCSA spectra. Comparison with another truncation method called mode truncation 
shows that the remarkably regular behaviour of the TCSA spectrum for large h in the case 
of the model (l2|), namely the presence of second flows, is not universal (i.e. not independent 
of the truncation scheme). The numerical calculations show that (l5|) provides a good 
approximation of the truncated spectra in both the TCSA and the mode truncation scheme. 
This is conflrmed by perturbative analytic calculations as well. The main difference between 
the mode truncated and TCSA spectra at large h seems to be explicable through the 
different behaviour of the function So{nc,h) in the two schemes. Difference between the 
So(nc, h) functions appears also in perturbation theory. We have shown analytically that in 
the mode truncation scheme the convergence of the truncated spectra to the exact spectra 
can be improved by (at least) one order in l/ric by the rescaling ([5]). This has also been 
shown in the TCSA scheme for low orders of perturbation theory in h. 

We have also given a quantum field theoretic discussion of the model (l2|) . In particular 
we have discussed (see Appendix [B|) the change of the boundary condition satisfied by the 
fermion fields as the coupling constant (or external boundary magnetic field) is increased. 
Such a change, which is emphasized in the literature, seems impossible naively — in our 
formulation at least. The paradox is resolved by the phenomenon that the fermion fields 
(more precisely their matrix elements between energy eigenstates) develop a discontinuity 
at the boundary if the coupling constant is nonzero. 

It is still an open problem to present an explanation of the validity of the approach ([5]). 
Within the framework of (l5l) the behaviour of the TCSA spectrum at large h, in particular 
the second fiow mentioned in the Introduction, is explained by the behaviour of sq and si 
at large h: Sq oc h, Si is bounded from above, therefore Si/sq tends to zero. It is a further 
problem to give an analytic derivation of this behaviour of Sq and Si. 

We shall investigate the scaling properties of Si/sq, Si and Sq and the TCSA and MT 
spectra in a further publication f33]. It is also interesting to extend the investigations to 
other perturbed boundary conformal minimal models, which show similar behaviour nu- 
merically to the model that we have studied. Certain results concerning other minimal 
models and scaling properties already exist [25], [20]; see also [18]. It is a further problem 
to classify the possible behaviours of truncated spectra at large h for various truncation 
schemes. Finally, the quantum field theoretic description of the model ^ could be devel- 
oped further and extended to case of massive particles. 
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A Distributions on closed line segments 

The necessary formulae for the Dirac delta 6{x) and step function Q{x) distributions on 
the closed interval [0, L] C R are the following: 



(98) 
(99) 
(100) 

(101) 

where / is a function defined on [0, L], and x e [0, L]. 

e{x-a) = if X <a , e{x - a) = 1 if x>a, (102) 



d{x 


— a)f{x) dx = 


/(«) 


if 


a e (0,L) 


6{x 


— a)f{x) dx = 


i/(0) 


if 


a = 


S{x 


— a)f{x) dx = 




if 


a — L 


5{x 


— a)f{x) dx — 





if 


a^[0,L] 



where a e R, 



d^e{x -L) = 26{x - L) (103) 
(9^0(,T -a) = 5{x -a) if a G (0, L) (104) 
d^Q{x-a) = Q if a^(0,L], (105) 



where x e [0, L] . 



^ exp[iA;(a; - x')] = 2L6{x - x') (106) 
exp[iA;(x + x')] = 2L[S{x + x') + S{x + x' - 2L)] , (107) 

where x, x' e [0, L]. 
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B Extension to Section [2 

In this appendix we present further details for Section [2l 

The anticommutators of the b{k) can be obtained in the following way: we define a 
scalar product on the classical complex valued solutions of the equations of motion: 

(^i,^2,a2|0i,02,&2) = / + V^2*02]dx + 2La;62 • (108) 

Jo 

This product should be calculated at a fixed time. Using the equations of motion it can 
be shown that the product is independent of this time. 

The essential properties of this scalar product are that it is defined by a local expression 
and that the n{k) are orthogonal with respect to it: 

{n{k,)\n{k,)) = 5fc,_fc„o(2L + . (109) 

Ki 

The creation/annihilation operators can be expressed in the following way: 

{nik)\i^^,^^I2,A2)) = bik){nik)\nik)) (110) 
Using the formula (11081) and the anticommutation relations (JG])-® we get 

m).ffe)} = 4,...o^^^^^^|^. (Ill) 

We remark that the above scalar product technique is also suitable for free fields on the 
half-line or in the usual full Minkowski space without boundaries in arbitrary spacetime 
dimensions. 

The nonzero matrix elements of the fields are 

(g|(v^,(x t),^2(x,t),A2(t))|P) ^ (113) 

^/\{P\P){Q\Q)\ 

where \Q) = h{k,)h{k2) . . . fe(^n)|0,), \P) = h{k^) . . . b{kjb{k)b{km+i) . • . &(^n)|0,), 

ALk 

n 

(QiQ)=n/(^^)- (115) 

i=l 
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^= E 7^^W^(-^)- (116) 

keS, k>Q 



The Hamiltonian operator can be written as 

k 

W) 

The following formula can be written for |0/i): 

N\Qn) = Jim e-°^|t;) = \{ - 77^^^)^-^)) 1^) , (117) 

k&S, k>0 \ ) / 

where is a normalization factor. The second equation on the right-hand side can be 
verified directly using the following formulae: [b{ki)b{—ki),b{k2)b{—k2)] = if /ci 7^ /c2 
(and ki,k2 > 0), and (6(A;)6(-A;))" = f {kY-%k)b{-k) . 

We remark that we have not given a mathematically completely rigorous proof that 
( l43l) satisfies (l6j)-(l9|), but we think that this would be possible. 

The expansion ( l43l) and (JED-Q imply that the following formulae hold: 



f{k)Mk){x, 0)M~k){y, 0) = 4L5{x - y) (118) 

kes 

J2fik)Mk){x,0)M-k){y,0)=4L6{x-y) (119) 

k£S 

J2 f{k)MkKx, 0)M-k){y, 0) = -4L[6{x + y) + 6{x + y - 2L)] (120) 

s 

J2 f{k)a2{kmM-k){x, 0) = (121) 

k£S 

J2 f{k)a2{kmM-k){x, 0) = (122) 
3 

J2fik)a2{k){0)a2{-k){0)=2. (123) 



kes 



fees 



kes 



These formulae are generalizations of fllOGp . (11071] . 
Using the formulae 



a{k) = -^^[ [ e''=^'^i(x,0) dx- f e-''=^^2(a;,0) dx] (124) 
2v2L Jo Jo 

Ai = ^[[ ^i{x, 0)dx- [ ^2(x, 0) dx] (125) 
^-^ Jo Jo 
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and (14311 we obtain the following relations: 

A-^E^M^^^ (127) 

. / N w, N — isin(A;L) v^w,n sm( kL) 

t^s t^s ^kLt^n{kL) 

Using (11101) we get the relation 



fc'Gf Z 



isin(A;L) , , , , , 

+ ^^^2(0) . (129) 

(I126l) - (ll29p can be regarded as Bogoliubov transformation formulae for this model. 
In the /i — > limit 

b{k{n,h)) -V2a{k{n,0)) n>l (130) 

b{k{0,h)) + b{k{0,h)y -Ai (131) 

i{b{k{0,h))^ -b{k{0,h))) ^ ^2(0) (132) 

\0h) ^ \v) . (133) 

The following boundary conditions are satisfied: 

(Ei|^i(0,t) + ^2(0,t)|^2) = {E,\^,mt) + ^2mt)\E2) = (134) 

and 

lim(Ei|^i(x,t) + ^2(a;,t)|^2) =0 (135) 

x—>0 

\im{Ei\d,^i{x,t) - d,^2{x,t)\E2) 

= 16Lh^{Ei\^2iL,t) - ^i{L,t)\E2) (136) 
lim lim(Ei|$i(x,t) - ^2(a;,t)|^2) = (137) 

h— >oo X— ►L 

lim lim(Ei|^i(x,t) + *2(a;,t)|^2) ^0 , (138) 

where and |i?2) are eigenstates of H. The boundary conditions ( 11341) are the same 
as those satisfied by the free fields. On the other hand, (I135p and (|136l) are similar to the 
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boundary conditions written down in [23l [271 I26]. From the point of view of the boundary 
conditions one can say that the perturbation Hj induces a flow from the boundary condition 
\imx^L{Ei\'^i{x,t) + '^2{x,t)\E2) = to the boundary condition \imx-^L{Ei\'^i{x,t) — 
^2(2;, t)|i?2) = and the boundary condition on the left-hand side remains constant, which 
is in accordance with the literature (see e.g. [23]). The boundary condition t) + 

'^2{L,t)\E2) = is called free spin boundary condition in the literature (e.g. [23j) and 
(£'i|\E'i(L, t) — '^2{L,t)\E2) = is called fixed spin boundary condition. 

We remark that from l\135\\ and (11361) and the bulk equations of motion (dt + dx)"^! = 0, 
{dt — dx)'if2 = equation f l38l ) can be recovered. 

We define the fields 

$i(x, t) = J2 -b{k)e"'^'-''^ ^2{x, t) = J2 KkW''^'^"^ . (139) 
kes kes 

They satisfy the following equations: 

(v&i - ^2)iL,t) = Y^^-('^2 - ^i)iL,t) (140) 
Mt) = -^{^i + ^2)mt) (141) 



Ho = -^ dx $i(a;, 0)9,<l>i(a;, 0) + — / dx $2(x, 0)a,$2(x, 0) - -hHi . (142) 

Note that the energies of the modes are not in -al^jL^ so one cannot conclude that the sum 
of the first two terms in p42l) equals to H. 

The above equations suggest how to describe the model discussed in Section [2] and in 
this appendix as a perturbation of the h ^ 00 limiting model. This description is discussed 
Section [3l The boundary conditions in this case are 



and 



(Ei|<l>i(0,t) + <f2(0,t)|E2) = {E^\^,iL,t) - $2(L,t)|E2) = (143) 



lim(Ei|<l>i(x,t) + <l>2(x,t)|E2) =0 (144) 
3;— >o 



lim{Ei\dx<l>i{x,t) - dx<l>2{x,t)\E2) 



x—^L 



16Lg 



(Ei|$i(L,t)-$2(L,t)|E2) (145) 



lim lim(Ei|$i(x,t) + $2(2;,t)|E2) = (146) 
lim lim(Ei|$i(x,t) - ^2{x,t)\E2) ^ , (147) 

g^oo x^L 



where \Ei) and [£'2) are eigenstates of if. 
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C Bethe-Yang equations 



The Bethe-Yang equations can be used to give a description of the spectrum of models in 
finite volume which have factorized scattering in their infinite volume limit. The Bethe- 
Yang equations for relativistic models defined on a cylinder are exposed, for example, in 
[HH |35l |36j ; for models defined on the strip they are written down in [371 [38]. It should be 
noted that the Bethe-Yang equations usually give approximate result only. 

In the case of the model that we study the ingredients of the Bethe-Yang description are 
the following: there is a single massless particle with fermionic statistics, the two-particle 
S-matrix is a constant scalar S{k) = —1, where k is the relative momentum. The refiection 
matrix on the left-hand side can be read from (11351) . it is Riik) = —1; the reflection matrix 
on the right-hand side can be read from (jl36p . it is 

16Lh'^ + ik , , 

The transfer matrices for A^-particle states are scalars: 

T^{ki,k2,...,k^) = RL{h)Rn{h) II S{h + k,) II S{h-k,) = -Rn{h) , 

i = l...N, (149) 

where fci > /c2 > • • • ^Af 7^ 0. This very simple form is the consequence of the simplicity of 
the S-matrix. The Bethe-Yang equations for the momenta ki, k2, . . . kN of the A^-particle 
states take the form 

e'^'^^^T^ik,, k2,..., k^) = e^^'^^il^j^ = 1 , t = l...N. (150) 

The total energy of an A^-particle state in the Bethe-Yang framework is E = ^i- (|150p 
can be rewritten as kiLta.n{kiL) = lQL?h?, i = 1 . . . N, which has the same form as ( i38l l. 
This means that the Bethe-Yang description reproduces the result of Section 12.21 for the 
spectrum exactly. 

The reverse model is similar, one can read from (11441) and (I145P that 

and the Bethe-Yang equations for the momenta can be written as 

kiLtan{kiL) = ^ , i = l...N, (152) 

which has the same form as (1621) . i.e. the result of Section [3] for the spectrum is reproduced 
exactly. 
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D Power series expansion of the energy levels 



The eigenvectors of Ho suitable for Rayleigh-Schrodinger perturbation theory are those 
introduced in (i2T]l . Degenerate perturbation theory has to be used. 
The nonzero matrix elements of Hj in Section 12.21 are the following: 

-{Qv\Hj\Qv) = {Qu\Hj\Qu) = 2 (153) 
{Qu\Hi\Pv) = {Pv\Hi\Qu) = 2y2(-l)"+"*(-l)'=^/" (154) 
{Qv\Hi\Pu) = {Pu\Hi\Qv) = -2V2(-l)"+'"(-l)^^/^ , (155) 

where P = a{ki)a{k2) ■ ■ ■ a{km)a{k)a{km+i) ■ ■ ■ a{kn), Q = a{ki)a{k2) . . . a{kn). 
We remark that certain perturbative calculations were also done in [39] . 



The eigenvalue of the state starting from a{^^)a{^^) . . . a{^^)\w) at /i = is 
(N, + N2 + --- + Nr)n 

\ ^ ' n=l / n=l 

where w stands for f or u, e = — 1 if = t>, e = 1 if w = u. rim = in the MT scheme 
is the truncation parameter introduced in Section SH]) , -{N^ + N2 + --- + Nr) 

in the TCSA scheme {ric is the conformal truncation parameter), and rim = C)0 in the 
non-truncated case. 

In the non-truncated case the coefficient of h"^ is ultraviolet divergent and should be 
regularized. The TCSA and the mode truncation both provide a regularization. 

We remark that the above formulae show that the truncated energy gaps converge to 
the non-truncated energy gaps as l/ric- 
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